

Iliill# 


V \ ; • ; 

Mi :: 


W-m3M§MUJ ■ : ill lllijliiii 

' ■ ■ i !1IS1 !!|Ii!I1!:I 




NASA TECHNICAL 







tlilSIl II : iliill 

■ liiiii ■ 

' : : ■ 1 


' : :;s sasssssii: 




by Joseph Norwood, Jr. 


WMmmSimlBaHmBSM. Va. 


CAUSES OF MAGNETIC-FIELD 
NONUNIFORMITY AND ITS EFFECT 
ONfcABIIi®l©F THETA-llICH PLASMA 




1 


gjli. 


I, 


NASA TR R-303 


CAUSES OF MAGNETIC-FIELD NONUNIFORMITY AND ITS 


EFFECT ON STABILITY OF THETA-PINCH PLASMA 


By Joseph Norwood, Jr 


Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 




CONTENTS 


Page 

SUMMARY . . . . . 1 

INTRODUCTION 1 

SYMBOLS 2 

INSTABILITIES OF THETA PINCH 6 

Drift Instability 6 

Rotational Instability .... ........................... . 11 

Other Instabilities 14 

SOURCES OF FIELD NONUNIFORMITY IN A THETA-PINCH COIL ....... . 18 

Collector Geometry 20 

Plasma Shape . 22 

End Effects ............... 23 

Coil-Collector Discontinuity ............................ 27 

Axial Current Distribution 29 

METHODS OF PLASMA STABILIZATION 30 

THE FIELD-SHAPING EXPERIMENT 34 

Description of Experimental Apparatus ...................... 34 

Use of Magnetic Probes ......................... 37 

RESULTS AND DISCUSSION 44 

Field-Shaping-Experiment Results . . . 44 

Stabilizing Effect of Flux Concentrator on Drift in Full-Scale Theta Pinch .... 55 

CONCLUSIONS 57 

REFERENCES 59 


iii 



CAUSES OF MAGNETIC-FIELD NONUNI FORMIT Y AND ITS 


EFFECT ON STABILITY OF THETA- PINCH PLASMA* 

By Joseph Norwood, Jr, 

Langley Research Center 

SUMMARY 

The various instabilities which beset the theta-pinch plasma are described in some 
detail, and the importance of the role of magnetic -field uniformity is demonstrated. 

Sources of the field nonuniformity are delineated and discussed. These include the 
collector-plate geometry, the nonlinear effect of the plasma, end effects in the coil, dis- 
continuity in gap spacing as the current flows into the coil, and effects of the alteration in 
axial current distribution due to skin effect. 

Next the various methods for stabilizing the plasma by eliminating field gradients 
are described, and a new method involving the use of a flux concentrator is proposed. The 
advantages of this method are compared with advantages of other techniques that have been 
used. Also, results of quantitative measurements of the field uniformity in an exact 
1/2-scale model of the Langley theta pinch are presented. These results include flux- 
concentrator and collector-plate-extension effects as well as the effect of the volume of 
the collector-plate structure which excludes high-frequency flux due to skin effect. 

Finally, experiments which suggest that full-scale theta pinches without extensions are 
capable of being stabilized by flux-concentrator techniques are described. 

INTRODUCTION 

The central problem in the confinement and heating of high-temperature plasmas is 
the production of a stable equilibrium, or, more realistically, the elimination or suppres- 
sion of those instabilities whose growth time is less than the attempted time of confine- 
ment. The theta-pinch geometry, in which a plasma column is confined and compressed 
by an axial magnetic field that is sinusoidal in time, has attracted considerable attention, 
and work on this configuration has been pursued at various laboratories throughout the 
world (ref. 1). The most dangerous of the instabilities which beset the theta-pinch plasma 

The information presented herein was included in a dissertation entitled "Stabiliza- 
tion of Plasma in a Theta Pinch" offered in partial fulfillment of the requirements for the 
degree of Doctor of Philosopy in Physics, Stevens Institute of Technology, Hoboken, 

New Jersey, 1967. 




have been shown to be caused by deviations from uniformity of the magnetic field (ref, 2) 
and by finite resistivity effects (ref. 3). The purpose of this paper is to point out the 
various causes of field nonuniformity in a theta pinch and to discuss the effect of such 
nonuniformity on the stability. Quantitative measurements of the uniformity of the vac- 
uum magnetic field in a 1/2-scale model of the Langley magnetic compression experiment 
(i.e., in a 1/2-scale model of the Langley theta pinch) are made for various conductor 
geometries. The effect of a flux concentrator on the field uniformity is compared with 
the effect of a collector -plate extension and the presence or absence of part of the flux- 
excluding volume of the collector-plate structure. Tests on the full-scale Langley theta 
pinch which demonstrate reduction of the drift velocity through the use of prototype flux 
concentrators are described. Probe measurements on a 1/5-scale model of the Naval 
Research Laboratory (NRL) theta pinch (ref. 4) were made for a restricted case with 
much less accuracy than was obtained in the present investigation. 

SYMBOLS 

A area of hole through which particle losses occur, meters2 

a radius of coil in conformal mapping, meters 

B magnetic flux density, teslas 

b width of transmission line, meters 

c speed of light, meter s/second 

c v specific heat per unit volume, joules/kilogram-meters^ 

d extension gap spacing, meters 

E electric -field strength, volts/meter 

e electronic charge (with power indicated by superscript u or w), coulombs 

F force, newtons 

g acceleration of plasma, meter s/second2; acceleration of gravity, 

9.8 meter s/second2 


2 



g(x,y) 

h 

I 

j 


Jz 


h 

k 

L 

l 

M 

m 


m. 


N 

n 

q 

R 

r 


r l 


r 2 

T 


spatial function of x and y 

gap spacing of transmission line, meters 

electrical current, amperes 

electrical current density, amperes/meter 2 

axial Hall current density, amperes/meter2 

azimuthal current density, amperes/meter2 

wave number, meter”! 

inductance, henries 

extension length, meters 

plasma mass per unit axial length, kilograms/meter 

mode number 

ion mass, kilograms 

number of particles 

particle density, meter “3 

charge, coulombs 

radius of curvature of field lines, meters 
radial distance from coil axis, meters 
inner radius of plasma column, meters 
outer radius of plasma column, meters 
temperature, degrees Kelvin 
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t 


time, seconds; thickness of coil wall, meters 


U = fi/k 

u real part of complex variable w 

v imaginary part of complex variable w 

v^ thermal ion velocity, meters/second 

V|| ion velocity parallel to B, meters/second 

v^ ion velocity perpendicular to B, meters/second 

Vg rotational velocity, meters/ second 

v mean ion velocity, meters/second 

v th mean thermal ion velocity, meters/second 

w complex variable; drift velocity, meters/second; extension width, meters 

W E zero-order drift, meters/second 

X radius of curvature of coil end, meters 

x real part of complex variable z 

y imaginary part of complex variable z 

Z coil axis; nuclear charge 

z complex variable; axial distance from the coil center, meters 

a probe discrimination coefficient 

(3 ratio of kinetic energy density to sum of magnetic and kinetic energy densities 

6 thickness of plasma boundary layer, meters; skin depth, meters 
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e 0 vacuum permittivity, farads/meter 

8 incremental angle, radians 

k Boltzmann constant, 1.3805 x 10~23 joule/degree Kelvin 

p m magnetic moment per unit length of plasma column, amperes/meter 

(i Q vacuum permeability, henries/meter 

p mass density, kilogram s/meter 3 

a electrical conductivity, reciprocal ohms/meter 

t plasma implosion time, seconds; acceleration time constant; duration of 

current pulse, seconds 

0 ion cyclotron frequency, second - 1 

co angular frequency of coil current 

cog reciprocal of instability growth time, second - * 

Subscripts: 

max maximum 

o compression field 

r radial component 

T trapped field 

z z- component 

r at time of plasma implosion 

00 at r = 0 and x = 00 

Dots over symbols indicate derivatives with respect to time; arrows over symbols 
indicate vectors. 
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INSTABILITIES OF THETA PINCH 
Drift Instability 

In many theta-pinch experiments the plasma column is observed to drift transversely 
in the plane of the collector plates either toward, as is usually the case, or away from the 
slot. The force F is due to interaction between the magnetic moment of the plasma 
column JI m , which is assumed to have only a z-component, and gradients of the confining 
magnetic field Bo, as indicated in the following equation: 

F = -(M m - V)§ 0 (1) 

The magnetic field in the coil (see fig. 1) is taken to be 



Figure 1- Magnetic-field model for drift calculation. 
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where B-p, the trapped magnetic field, may be either positive or negative. From equa- 

tion (1), F r = fj. m however, since the field in the coil for r > r 2 is due only to 

currents in the external conductors, the curl of B 0 vanishes and 8B o ,r/ 0 z maybe 
replaced by 9B 0 ,z jdv. The transverse acceleration of the plasma column may then be 
written as 


r 


dB, 


Mm 


o,z 


dr 


M 


where pt m is the magnetic moment per unit length of the column and M is the plasma 
mass per unit axial length. The magnetic moment per unit length is given by 


Mm 



(3) 


Using Ampere's law for the current density 


and equation (2), one finds 


1 dB z 
M 0 dr 



Thus the transverse acceleration of the plasma column becomes 


7r ^o ~ ^T ^o,z r 2** ~ r l^ 
3 ix Q M dr r 2 - r j 


(r < r 2 ) 


(4) 


(5) 


With 

1 dB o,z _ _ 1_ 

B 0 dr R 

where R is the radius of curvature of the field lines and where B 0 = B max sin o>(t + t) 
^with B max being the peak magnetic flux density, w the angular frequency of the coil 
current, and r the plasma implosion time), equation (5) integrates to 
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2 3 3 / r 

_^max 2 l_[Xp_(cos 2w(t + t) - cos 4 cut] +— (t - T)(t + 3 t) 

3m 0 MR r 2 - r 1 \2co\4co L J 2 

+ - - - sin 4o;t\ + — — f ifsin cu(t + r) - sin 2o)rl - (t - t)cos 2cot 

2 J wB max ^col J 


( 6 ) 


This equation is derived with the assumption that the linear mass density M and 
radius of curvature of the field lines R are constant. Considering the plasma implosion 
time t to be negligibly small and letting be zero with r ^ consequently going to 
zero, one finds 


B 2 
D max 

1( 

6p 0 nmjRu> 

4w 


cot 2 

2 


(7) 


By taking values of B max , n, m|, and w appropriate to the Langley theta pinch and 
using a value for R measured in the Naval Research Laboratory experiment (ref. 4) , 
which is geometrically similar, the curve shown in figure 2 is plotted from equation (7). 
The constancy of R and n is confirmed by plotting the measured displacement of the 

plasma column as a function of ^-(cos 2o>t - 1) + and comparing the results with a 

straight line. This constancy is verified by the moderate point scatter shown in figure 2. 
The measured points were taken by streak photography with a rotating-mirror camera 
looking at the plasma in the median plane of the coil from the bottom — that is, recording 
plasma motion in the horizontal plane. The streak photograph is shown in figure 3. This 
shot was made with a main-capacitor -bank energy of 335 kilojoules and a filling pressure 
of 100 millitorr (13.3 N/m2) of hydrogen. No bias magnetic field was employed, but 
because bias fields are often used, the trapped magnetic field Bx representing a bias 
field is included in the theory. 

The agreement between theory and experiment is not as good at higher energies, 
as can be seen in figure 4, which was plotted from measurements of the streak photograph 
shown in figure 5. This shot was made with a filling pressure of 13.3 N/m2 of hydrogen, 
no bias field, and an energy of the main capacitor bank of 1 megajoule. The lesser agree- 
ment stems from the fact that the particle density n and the radius of curvature of the 
field lines R were assumed to be constant. The assumption of a constant particle den- 
sity is certainly poor at high energies since, without magnetic mirrors or a reversed 
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bias field, the plasma rapidly escapes from the coil ends and the particle density at the 
axial center can be expected to decay markedly while the plasma column drifts. 

The effect of such end losses can be included in the theory. Equations of state for 
an adiabatigally compressed plasma with end losses have been derived in reference 5, 
and the following equation obtained from that reference describes a j3 = 1 plasma whose 
ions stream out the ends with mean thermal velocity: 


T = 




( 8 ) 


In this equation, B = B max sin wt and B T = B max sin cur. The particle loss rate, which 
is just the efflux rate through a hole of area A (ref. 6) , is represented by 


dN _ o nv th 
dt 4 


(9) 


where 


and where the factor 2 is included because losses occur at both ends of the tube. Equa- 
tions (8), (9), and (10) constitute a set which can be solved numerically for n. The result 
can be used in equation (7) to yield the loss-corrected drift. This process has not been 
carried out, however, since the necessary input parameters (temperature, volume, and 
especially hole area) are not known to an accuracy necessary to justify this degree of 
elaboration in the theory. 

It should be noted that the drift of a theta-pinch plasma due to field gradients is not 
truly an instability of some equilibrium state but rather a nonequilibrium situation, the 
correction of which depends on the reduction of the field gradients to the point where the 
drift times are long compared with the half -cycle period of the compression field. 

Rotational Instability 

Several mechanisms leading to rotation of a magnetically confined plasma have been 
postulated. The instability by which this rotation is observed in a theta pinch takes the 
form of a fission of the plasma column into two strands which separate, go to the walls, 
and thus destroy the confinement. This effect is depicted in figure 6 by a line drawing 
taken from reference 7. The rotational instability has been observed photographically 
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1 

Figure 6.- Line drawing of streak showing m = 2 
rotational instability. 


both with side -on and end-on streak photog- 
raphy (refs. 8 to 11) and with end-on framing 
photography (refs. 8 and 11). 

The effect is usually interpreted as an 
m = 2 Rayleigh-Taylor instability caused by 
the centripetal acceleration of the plasma 
column in the presence of the confining mag- 
netic field. The single-fluid hydromagnetic 
description of this instability (ref. 12) leads 
to the dispersion relation 

2 = -kg (11) 


where cog is the reciprocal of the growth time of the instability, k is the wave num- 
ber, and g is the acceleration of the plasma. This relation predicts unstable solutions 
for all wave numbers. With the inclusion of finite ion cyclotron frequency and radius 
(see refs. 13 to 15) in the theory, the following equation for the stability criterion (given 
in the notation of the present paper) was obtained in reference 16: 



( 12 ) 


The thickness of the plasma boundary layer 6 is defined as — ^ — . Equation (12) is 

dn/dr 

ni/'p 

valid only for jS = 0, where j3 = — — — . However it can be shown by using the equations 

B 2 /2p 0 


of reference 16 that the correction for finite /3 is negligible. The inequality may be 
solved for plasma acceleration to yield 





J 


where U = O/k. For the rotational instability, k ~ raj Y2 where m is observed to be 
2. By taking B ® 5 teslas and considering a hydrogen plasma with v = 500 km/sec 
and 6 = 1 centimeter, the range of unstable angular velocities is found to be 
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5 x 10 6 = co i 3. 9 x 10 8 


(14) 


which is in agreement with the observed unstable angular velocities (refs. 8, 9, 11, 17, 

18, and 19). 

The mechanism by which the plasma acquires angular momentum is not altogether 
clear. A large number of feasible mechanisms have been proposed to explain the phe- 
nomenon. These theories are discussed comparatively in reference 20. They can be 
classified according to how the reaction torque is taken up. Purely electromagnetic 
angular momentum as the reaction absorber (ref. 21) can be shown to be too small by 
four to eight orders of magnitude in a large theta pinch although it may explain rotation 
in a low-density mirror machine. Roberts and Taylor (see refs. 8 and 22) have suggested 
that the reaction is taken up by circulating currents at the wall and that the plasma is 
driven essentially by the electric field in the sheath. In reference 20 Haines has modified 
an earlier suggestion by Jensen and Voorhies (ref. 23) to show that viscous momentum 
transfer to the walls by ion- wall collisions may account for the rotation. Bostick and 
Wells (ref. 24) have suggested that plasma is expelled as vortex rings from either end 
while the remaining plasma, the toroidal velocity of which is opposite to that in the rings, 
coalesces in the center. Velikhov (ref. 25) has suggested a radial division of the plasma 
into oppositely rotating layers. The best supported theory is that of Haines in refer- 
ence 26, where a model in which the reaction is absorbed by external conductors via the 
magnetic stresses is proposed. This model depends on the presence of a radial compo- 
nent of magnetic field with a nonvanishing azimuthal gradient. The ideas can be summed 
up in a simplified analysis written by the present author. An axial Hall current density 
j z is generated by the interaction of the azimuthal current density j q with the radial 
component of the magnetic field. The Hall current then interacts with the radial magnetic 
field to produce an azimuthal force on the plasma column proportional to B r ^. A simple 
statement of the generalized Ohm’s law and Newton's second law yields 


and 


ie B i 

ne 


v 0 B r 


dv fl 

P W-^ B r 


(15) 


(16) 


The first term on the right-hand side of equation (15) describes the axial Hall electric 
field. Equations (15) and (16) are combined to yield 
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(17) 


dv 0 

dt 


oB 2 

+ 

P 


V 0 = 


°ie B r 2 

nep 


which, with vg = 0 for t = 0, gives for the rotational velocity 


v e 




(18) 


Thus the limiting rotational velocity is appropriate to the velocity associated with diamag- 
netic current when j z vanishes. If jg is reduced to less than nevg, the Hall current 
flows in the opposite direction and decelerates the rotation. Thonemann and Kolb (see 
refs. 3 and 20) have analyzed, with a computer, the full two-fluid equations including field 
distortion by the Hall currents and have found a time constant for the acceleration process 
given by 


T = -P-lJm ) 2 

aB r 2 Vne/ 


(19) 


A summary of the main features of these rotation mechanisms is taken from refer- 
ence 20 and shown in table I, with the notation changed to that of the present paper. 

Strong experimental evidence has been found to support the Haines model of rotation that 
involves the external conductors. It has been shown in reference 27 that an externally 
applied transverse magnetic field of only 0.015 tesla is sufficient to induce rotational 
instability. Thus, the rotational instability and the drift of the plasma column have both 
been shown to depend critically on the nonuniformity of the magnetic compression field. 

Other Instabilities 

The instabilities discussed in the two preceding sections constitute the most dan- 
gerous of the instabilities which beset the theta pinch in the sense that they limit the con- 
finement time to a quarter or less of the coil-current oscillation period. The theta pinch 
also exhibits a rich variety of less drastic instability phenomena which are now briefly 
discussed for the sake of completeness. 

A "wobble" instability was found in the Culham Laboratory theta-pinch experiment 
(ref. 28). This instability manifests itself as a spiral precession of the center of the 
plasma column. It has been suggested that wobble instability is due to slow rotation of 
the plasma column (below the threshold for the m = 2 instability) and to magnetic line- 
tying between the column and currents at the vacuum -chamber wall. Line-tying gives 
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an effective friction with the wall. As the plasma column drifts off the magnetic axis, 
some of the column spin becomes transformed into orbital angular momentum of the 
center of mass. The column is not observed to strike the wall within the confinement 
time. 

In reference 29 it has been predicted that a theta-pinch plasma with reversed 
trapped field would be subject to a resistive instability resulting in a "tearing"-mode 
breakup (breakup along current-flow lines) of the column into short rings. This process 
has been observed in reference 30 and is illustrated in figure 7. In experiments in which 
short coils are used, the current layer often collapses into a single ring (ref. 8), presum- 
ably because the tearing mode which grows fastest for long wavelengths does not have 
room or time to establish itself. For longer coils, as many as seven rings appear 
(ref. 30) and subsequently coalesce as axial contraction of the plasma column occurs. It 
has been suggested that vortex rings of the Alfven wave type (ref. 3 1) similar to those 
observed in references 32 to 36 may be present in the Langley theta pinch. Transverse 
striations are seen in the streak photographs of figures 3 and 5 (taken side-on). The 
initial striations are dependent on the magnetic field strength and the linear mass density 
of the plasma (ref. 37) and are strongly damped as a result of plasma heating by magnetic 
pumping (ref. 38). A series of striations of different period follow the initial striations 
and show much less damping. A microdensitometer trace of a streak photograph illus- 
trating these effects is shown in figure 8, where the film density is recorded along the 
center line of the plasma image. The plasma from which these data were taken had no 
trapped field; hence one would not have expected the tearing-mode instability to develop. 

It is unlikely that the observed series of striations of different period are associated with 
radial oscillations as are the initial oscillations, since no modulation of the plasma radius 
is evident. Rather, they are tentatively interpreted as magnetovortices generated by a 
resistive "rippling" mode (ref. 29) that propagate along the plasma column with Alfven 
speed. This interpretation has not been verified since poloidal flows necessary for con- 
firmation would have to be detected by probe measurements and such measurements are 
impractical in such a high-energy plasma (refs. 39 and 40). Velocity measurements, 
however, may be possible. 

In short coils, axial magnetic-line tension tends to concentrate the plasma into a 
single short ring which subsequently seeks to flip over so as to aline its trapped field with 
the magnetic compression field (see ref. 41). This instability has been analyzed in ref- 
erence 42 where plasmas with sufficiently high length-radius ratios are shown to be stable 
with respect to flipping. 

Azimuthal and axial striations have been observed in the photographs presented in 
references 18, 43, and 44. Also, Rayleigh-Taylor instabilities (flutes) are evident in end- 
on framing photographs in these references. It has been suggested that these long 
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serpentine flutes may be due to helical kinks in the axially alined striations. An analogy 
to a rope twisted under tension might be made; kinks appear as the tension is relaxed. 

SOURCES OF FIELD NONUNIFORMITY IN A THETA-PINCH COIL 

As the foregoing discussion points out, the most harmful of the plasma instabilities 
which occur in a theta pinch — that is, the drift and rotation — are caused and/or aggra- 
vated by nonuniformity in the magnetic compression field. There can also be effects due 
to lower electrical conductivity of the plasma caused by the presence of impurities of 
high nuclear charge. In reference 27 it has been shown that the presence of as little as 
0.1 percent of contaminant of high nuclear charge can cause explosive rotational instabil- 
ity. If the contamination is appreciable, the electron temperature (and hence the electri- 
cal conductivity which varies as T 3 / 2 ) will be seriously degraded by radiation. For high 
densities such that the ions are effectively coupled to cold electrons, the finite-gyro- 
radius stabilization of a rotating flute would be less effective (see ref. 15). In addition, 
the magnitudes of the torques which drive the rotation depend on the electrical conductivity 
(ref. 3). In the drift the effect of contamination is not so drastic. The magnetic moment 
can be altered somewhat by changes in the skin depth of the plasma and there can also 
exist secondary effects, which are discussed subsequently. The problem of minimizing 
contaminants of high nuclear charge is essentially one of vacuum technology and proper 
gas injection and is not discussed herein. Rather the present discussion is concerned 
with the causes of field nonuniformity in a theta-pinch coil. 

The requirement for connecting a large number of coaxial cables (up to 1000) to a 
single-turn coil about 1 meter in length dictates the geometry to a large extent. In the 
Langley magnetic compression experiment as well as in the NRL experiment, a pair of 
large (4.8 by 1.2 meters) current-collector plates are used with the cables connected to 
one of the long sides and the coil symmetrically attached to the other. This system is 
shown schematically in figure 9. Figure 10 is a photograph of the Langley magnetic com- 
pression experiment. The lead bricks on the collector plates are necessary to limit the 
amplitude of the bouncing of the plates when the current pulse passes through them. The 
plates are further restrained by 128 steel bolts. The coil is shown in the center of the 
photograph and in more detail in figure 11. Also shown in figure 10 are the rotating - 
mirror streak camera and a spectrograph. By considering this general configuration to 
be typical of large theta pinches, several major sources of magnetic -field inhomogeneity 
can be delineated. 
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Collector Geometry 


Figure 12 shows a resistive analog plot of the field lines to be expected in the col- 
lector plates of the Langley magnetic compression experiment. The two plates have been 
unfolded about the coil in order to create a planar surface. The equipotential pattern 
shown corresponds to magnetic field lines in the limit of dc operation, since the mag- 
netic field can be described in current-free space by a scalar potential which obeys the 
Laplace equation. The reproduction of the ac conditions occurring in the experiment 
(skin depth much less than conductor thickness) is technically unfeasible in the analog 
plot and hence the region near to and including the coil is not altogether accurately 
represented. 

The top plate is slotted perpendicular to the long sides, as shown in figure 12, in 
order to provide a high inductance path to inhibit destructively high currents in case of a 
short circuit between the plates. The holes are provided to bolt the plates together. 

The field pattern on the left side of the center line is obtained with current supplied 
only to the center third of the plates. The pattern on the right side is obtained with cur- 
rent supplied to the entire length of the plates. The field lines near the coil-collector 
interface show considerable curvature under both conditions — that is, the current has a 
nonnegligible component parallel to the coil axis. It is this "funneling" effect of the cur- 
rent and the consequent gradient of the field in the coil that probably contribute most to 
drift instability in large theta pinches. 

As pointed out in reference 2, the curvature may be aggravated by too large a coil 
feed slot. Since curved field lines in the coil may bend into it, this slot should be mini- 
mized in order to minimize the curvature. 

Another source of field gradient due to the collector geometry may be the large vol- 
ume external to the coil from which the field is excluded. The collector plates and their 
supporting structure together with the lead bricks, all on one side of the coil (see fig. 10), 
exclude the high-frequency field from this volume and thus one might expect some degree 
of azimuthal asymmetry of the field also in the coil interior. It is suggested in refer- 
ence 45 that the effect might be fairly large in the end regions. For long coils, as is 
shown subsequently, the ratio of magnetic field outside the coil to magnetic field inside 
the coil is on the order of only a few percent; hence the field stresses, which vary as B^, 
are almost four orders of magnitude less on the outside than on the inside. Thus devia- 
tions from homogeneity of the inside field far from the end (e.g., in the median plane) due 
to external azimuthal asymmetry should be less than 0.1 percent. In reference 46 the 
introduction of large metal plates at the ends of the coil to partially balance the azimuthal 
symmetry was noticed to cause some difference in drift rate. A similar configuration 
had been tried previously and no such effect observed, so the question of the influence of 
end plates on drift remains an open one. 
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Current supplied to 
center section only 


Current supplied to 


L-68-10,076 

Resistive analog plot of magnetic field lines in the collector plates of 
the Langley magnetic compression experiment. 












Plasma Shape 


When the plasma contracts axially such as to concentrate itself near the median 
plane of the coil, that portion of the coil has lower inductance and the driving current 
converges there. The effect is similar to the current funneling just discussed. The 
field lines have a curvature that is convex toward the coil slot, and the plasma may be 
expected to drift in that direction. 

The axial contraction is usually observed to be a strong effect only in plasmas 
employing a reversed bias field. As shown previously herein, resistive instabilities 
and/or vortex coalescence may serve to cause axial contraction. It can be shown that a 
tubular plasma corresponds to an equilibrium state and that tension of the axial field lines 
should not cause contraction. Thus an instability of some sort must provide a mechanism 
to trigger the contraction. The problem has been investigated in the Langley theta pinch 
by examining the plasma with the use of light pipes (fiber bundles) spaced every 15 centi- 
meters along the coil axis. The ends of the light pipes away from the axis form a con- 
tinuous slit which is focused on the photocathode of an image-converter streak camera. 
The monitoring system, which is similar to one used at Langley Research Center to study 
spark-gap jitter (ref. 47), is shown in figure 13. Figure 14 is a typical axial-light- 
intensity streak record. The sinusoidal shape is due to perturbation of the electron optics 
in the image-converter tube caused by the magnetic field of the discharge and actually 
allows direct comparison of the axially resolved light intensity with the discharge cur- 
rent. A trace shows that axial contraction begins roughly 5 ps after the start of the 



Figure 13.- Axial-light-intensity monitoring system. 
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L-68-10,077 

Figure 14.- Axial-light-intensity streak record. 


discharge. The ends appear to move in with a speed of approximately 5 cm/pts until the 
column has contracted axially to a length of approximately 40 centimeters. Illumination 
in the end regions due to end losses builds up until, near the end of the first half -cycle, 
the plasma streams out at the coil ends, follows the flaring magnetic field, and strikes 
the walls. Impurity line radiation greatly intensifies at that time. With large end losses, 
there is difficulty in determining whether it is axial contraction or an axially propagating 
rarefaction wave that is observed. The effect on the magnetic field is, however, the same 
The second half-cycle is largely overexposed as a result of the intense line radiation. 
These phenomena correspond rather closely with those observed in reference 27 for the 
NRL experiment. The plasma contraction is apparently limited to a minimum length of 
about a third or a quarter of the initial column length. The observed plasma radius does 
not seem to change in this time; thus for the high compression ratios employed (approxi- 
mately 100:1), the change in inductance should be moderate. In a short coil where the 
axial contraction is a stronger effect, the perturbation on the field may be significant. 

End Effects 

Because of the finite length of the coil, the field lines diverge and thus produce an 
azimuthally symmetric radial gradient which increases as one approaches the end of the 
coil from the center. This field pattern is not the same as for the dc solenoid where 
flux can permeate the coil wall. For the ac conditions (skin depth much less than radial 
coil thickness), the field pattern is approximately that depicted in figure 15, where the 
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constant-v lines represent the magnetic 
field lines. The coil is considered to be 
semi-infinite. A Schwarz -Chris toff el 
transformation yields the mapping from 
this two-dimensional geometry to the 
Cartesian coordinate system. The equa- 
tion for the Maxwell curves is found to be 
(ref. 48) 

z = + 1 + ewj (20) 

where 


Figure 15.- Maxwell curves obtained from equation 
z = f.(w + 1 + ew). 


with i = fl. Equation (20) can thus be written as 


z = x + ly 


w = u + IV 


( 21 ) 


x 


^1 + u + e u cos vj 


( 22 ) 


or 


y = + eUsin vj 


(23) 


These equations may be inverted by numerical iteration in order to find 8v/8y = g(x,y) 
from which B(r ,x) may be obtained. The results are shown in figure 16 where, with x 
being the distance into the coil from the end, the ratio of B(r = 0,x = 0) to 
B(r = 0,x = °°) - that is, B/Boo - is plotted as a function of r/a where a is the coil 
radius. Further calculations show that B/B^ approaches unity to within 0.1 percent at 
a distance of two radii from the end. This result supports the use of a semi-infinite coil 
length in the calculation instead of a finite coil length. Thus stabilization against drift due 
to the "minimum B” magnetic well is confined to the very ends of the coil and cannot be 
expected to be effective in large experiments where the coil aspect ratio (ratio of coil 
length to radius) is greater than 20. The fact that the calculated mirror ratio 
B(r = 0,x = 0) 

— — — — — agrees within a few percent with the measured value gives some degree of 

Boo 

confidence in the two-dimensional model. 
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Figure 16.- Magnetic contours at the end of a semi-infinite coil. 


As shown in reference 49, the plasma drift may be calculated from the point of view 
of orbit theory as well as by use of a magnetohydrodynamic model. The orbit-theory 
analysis in reference 49, although correct, involves currents and thus is somewhat a 
hybrid. The following analysis by the present author is therefore preferred. The equa- 
tion for the velocity of particle drifts that cause charge separation (obtained with the use 
of ref. 16) is 

w = 


where the time rate of change of the electric-field strength is given by 




(25) 


(26) 
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iUlll 

where 1 +— is the plasma dielectric constant, which is large compared with unity. 


e o B 


Thus equation (26) reduces to 


£ 



(27) 


The zero-order electric drift which acts on the plasma as a whole is given in reference 16 
by 


W E _ E X B 
B2 

Thus the following equation can be written: 

v i 2 W 

2 j B 



(28) 


(29) 


The term involving dB/dt has been neglected in equation (29) inasmuch as the reciprocal 
of the time for polarization of the plasma v fa Jr 2, where v^ is the thermal velocity of 
the ions and r 2 is the radius of the plasma column, is typically two orders of magnitude 
larger than the angular frequency of the magnetic flux density o>. With 


VB _ 1 
B R 


(30) 



and, for the assumption of a /3 = 1 plasma, 

■r2 

2n/cT = — 

2/i 0 


equation (29) can be written as 

» _ B 2 TrB 2 r 2 2 
2nm i p 0 R 2Mp 0 R 


which is, within a factor 1.5, the same as equation (5) with B r = r^ = 0. 


(31) 


(32) 


( 33 ) 
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The point of this analysis is that there exists a polarization in the drifting plasma, 
at a right angle both to the direction of drift and to the magnetic field, whose sign is 
determined by the direction of the curvature. In reference 50 it has been found that in a 
theta pinch of moderate length, considerable curvature can be tolerated without drift. 

This result has been verified in reference 51. The explanation is that gradients reverse 
in the end regions, particularly for those coils which feature magnetic mirrors. Thus the 
polarity of the charge separation in the end regions is opposite to that near the coil cen- 
ter. The electron conduction along magnetic field lines allows the fields to be short- 
circuited. (See fig. 17.) This effect has been demonstrated very dramatically in refer- 
ence 50 by isolating the end regions from the main body of the plasma with the use of 
glass plates. The plasma then drifts rapidly to the wall. 



(a) Top view, with polarization normal to the page. 


, jl . 4- 4- -f - 4- -f. j 



(b) Side view showing reversal of polarization due to inflection in 
curvature of field. 

Figure 17.- Polarization effect on drift. 


Coil -Collector Discontinuity 

In the high-frequency limit (i.e., with the skin depth 6 much less than the thickness 
of the conductors), the partition of the current between the inside and outside surfaces of 
the plates in a parallel-plate transmission line is a function of the gap spacing h and 
the width of the transmission line b. As pointed out in reference 2, a variation of the 
ratio h/b along the length of the transmission line causes current to spill around the 
ends in order to adjust the current ratio (ratio of current outside to current inside) to the 
value necessary for maintaining a constant skin depth. This effect is of interest herein 
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since the collector plates, coil slot, and coil form a transmission line with abrupt changes 
in h/b. 

As the current flows from the coil slot (h/b ~ 0.001) into the coil (h/b ~ 0.1), the 
value of h increases from 1 millimeter (at the slot) to 10 centimeters (the coil inner 
diameter). Part of the current flowing on the inner surface must then spill around the 
ends of the coil. 

The ratio of current outside to current inside may be computed from equations (22) 
and (23). The quantity 9v/9y at y = a (for v = it) is required (see fig. 16). Setting 
v = tt ± 6 , where 9 is small, and expanding to first order yields 

x=2(i +u -e u ) (34) 

and 

y =|ji ± ^(e u - !)] (35) 

for which 

_90_9v _ 7r/a 

9y 0y e u _ 1 

V=7T 

Equations (34) and (36) serve to determine the ratio of magnetic -field strength at any 
axial location in the coil. In order to find the ratio of total current outside to total cur- 
rent inside, one must integrate equation (36) with respect to x by means of equation (34). 
However since x = 0 corresponds to u = 0, equation (36) is singular at the origin. 
Physical arguments are advanced in the next section to determine the lower limit on this 
integral. As an indication of the magnitude of the current partition, the ratio of magnetic - 
field strengths can be compared at the median plane of the transmission line by setting 
x = b/2 in equation (34). The results, with values appropriate to the Langley experiment, 
are shown in figure 18. The effect on the field gradient in the coil due to this spilling of 
current around the ends is opposite to that induced by the collector geometry — that is, 
the gradient is directed toward the slot instead of away from it. 


(36) 
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Figure 18.- Divergence of the current from the inside to the outside of a parallel-plate 
transmission line for x = b/z. 


Axial Current Distribution 

The axial distribution of current in a high-frequency coil is not uniform as in the 
dc case. Equations (34) and (36) may be used to calculate the axial distribution. The 
current-density variation near the 
end of the coil is given in figure 19. 2 

In practice, the limiting value of the \ 

current density is determined by \ T otal i 

the coil material through the 

requirement that the heat generated 10 \ \ 

at the conductor surface by joule \\ 

losses not melt the coil material. arbitrary units 

This requirement is equivalent to io° - — 

specifying a minimum radius of \ 

curvature of the edge at the end of \ 

the coil. As shown in reference 52, 10 -i_ \ 

the limiting value of a sinusoidally \ 

varying, high-frequency magnetic \ 

field that can be tolerated at the \ . 

surface of a conductor can be given 10-5 10 ' 4 10 ‘ 3 io - 2 io - 1 io° io 1 io 2 

by the relation (in notation of pres- x 

ent paper) 


AT =• 


2/i 0 c v 


Figure 19.- Current-density variation near the end of the coil. 
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where AT is the temperature gain per half-cycle of the current oscillation, B max is 
the peak magnetic flux density at the coil surface, and c v is the specific heat per unit 
volume. For the beryllium -copper coil used in the Langley experiment the maximum 
tolerable field, with 1 half-cycle of operation assumed, is 54 teslas. In order for the 
Langley coil to be used with the full capability of the capacitor bank (12.6 megam peres), 
the edge of the coil end must be rounded to approximately 0.5 millimeter. If the edge is 
not rounded, the field itself will round the corner to conform to the 54 -tesla field contour. 
This effect has been recorded photographically. During about the first hundred shots of 
the coil, sparking from the edges was observed. As the coil edges eroded, the sparking 
diminished until there was no point on the coil where the surface field exceeded the criti- 
cal value. 

Thus the partitioning of the total current between the inside and outside surfaces 
depends not only on h/b but also on the value of the total current. The determination 
of the current ratio involves numerical solution of an integral of the form 



where t is the thickness of the coil wall, 6 is the skin depth 


X is the 




radius of curvature of the coil end ^itself a function of ^total)’ anc * ^ * s the constant to 

be determined. The integration must be performed for u > 0 (the outside surface) and 
u < 0 (the inside surface). 


As the current flows from the collector through the feed slot into the coil, an 
increasing percentage of the total current is channeled away from the median plane, 
inside as well as outside. This facet of the current redistribution phenomena reinforces 
the effect of current transfer described in the preceding section. The field gradient is 
directed toward the feed slot in opposition to the effect of current tunneling. 


METHODS OF PLASMA STABILIZATION 


Evidence that major instabilities of theta-pinch plasma are caused by nonuniformity 
of the magnetic field as well as by the presence of impurities has been presented. The 
problem of reducing impurities involves advanced vacuum technology and the design of a 
high-purity-gas filling system such as the palladium or nickel leak which can be used to 
admit hydrogen or deuterium. This aspect of the problem is, however, not within the 
scope of the present work. The source of the field nonuniformity has been discussed. 

The major contribution to drift toward the feed slot (which is the dominant direction of 
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drift) seems to be the current funneling caused by the collector geometry. Axial con- 
traction of the plasma column probably contributes also. The magnitude of the plasma- 
axial-contraction effect is difficult to ascertain since there are also end losses and the 
situation has not yet been described accurately by theory. Drift away from the feed slot, 
which is observed to be a smaller effect 
(ref. 2), is normally caused by current 
redistribution. Drift away from the feed 
slot has been observed in the Langley 
experiment in a plasma with reversed 
trapped field. The drift is reversed before 
the outer wall is reached and motion toward 
the feed slot ensues. (See fig. 20.) This 
effect seems to occur only when a reversed 
trapped field is used. The most likely 
explanation is that instabilities connected 
with the reversed-field configuration some- 
how enhance the motion of plasma toward the ends where it apparently is strongly con- 
centrated, with an effect just opposite to that observed when the plasma contracts axially. 
The turn-around is caused by the expulsion of the plasma at the ends and/or by forces 
that arise as the plasma approaches the conducting wall of the coil (see refs. 53 and 54). 

It has been found in the Langley experiment that the rotational instability does not 
prove troublesome during the first half-cycle. After the plasma strikes the wall as a 
result of drift or the oscillation of the field, rotation is noted immediately on the following 
half-cycles and the plasma is violently unstable. Only the first half-cycle is used for 
observation; thus the drift constitutes the only disturbance which seriously limits the 
experiment. 

The practical approaches to the solution of the instability problem in theta pinches 
are typified by the methods employed at NRL and at Culham Laboratory. Initially, NRL 
(see refs. 4 and 55) and Los Alamos Scientific Laboratory (see ref. 56) both used an Ioffe 
(ref. 57) hexapole magnetic field in an attempt to center the plasma column. The radial 
magnetic field associated with the Ioffe bars strongly excited the rotational instability, 
and thus this method had to be abandoned. Next, NRL built an extension to its collector- 
plate system, with a width equal to the coil length and with a length of 1 meter. This 
extension was found to reduce the measured value of AB z /b z of 0.1 percent by approx- 
imately a factor of 10 (see ref. 27). There are strong disadvantages to such an extension, 
however. 

The coil inductance in a large theta pinch is typically about 7 nanohenries. In order 
to get a reasonable fraction of the bank voltage applied to the coil, the parasitic inductance 



L-68-10,078 

Figure 20.- Streak photograph showing retrograde drift in 
plasma with reversed trapped field. 
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must be kept down to about 2 nanohenries. The parasitic inductance added by the 
collector-plate extension is 

T M d 

w 


( 39 ) 


where l is the extension length, w is the extension width, and d is the extension gap 
spacing. For a square extension ( l = w) and a 1-millimeter gap spacing, the inductance 
would be about 1 nanohenry. Thus the voltage at the coil (upon which the ion heating by 
shocks depends) would be seriously degraded. In addition, the channeling of large cur- 
rents would produce extremely high mechanical forces. The force between the plates 
would be 


r M 0 I h 

2w 


(40) 


With a mass loading M on the plates (usually consisting of lead bricks), the amplitude 
of the bouncing - that is, the maximum gap spacing — can be computed to be 


■hnax 


MoI 2 Zt2/p 0 I 2 ^ 
4wM \2wMg 


+ 1 


(41) 


where g is the acceleration of gravity and r is the duration of the current pulse 
(assumed to be a square pulse). For a 1-meter-square extension with 10 megamperes 
of current for 20 jbi s and a mass loading of 10 000 kilograms, the amplitude of bouncing 
is approximately 1 millimeter. This bouncing causes severe hammering of the insulation. 

It should be noted that although the current-funneling effect is adequately compen- 
sated, a simple extension does not cancel field nonuniformity caused by the coil-collector 
spacing discontinuity. In this respect the approach at Culham Laboratory (see ref. 2) is 
somewhat more sophisticated. A short extension which houses an inductive lens is used. 
(See fig. 21.) The inductive lens is an interchangeable section of the collector which has 
an increased inductance. The inductance varies across the width of the tab extension in 
such a way as to compensate for the convergence of the current paths in the coil. The 
objections which have been raised for the simple extension apply to this method as well, 
but the advantage of this technique lies in its flexability. 

The solution being pursued in the present paper involves the use of flux concentra- 
tors (ref. 58). A flux concentrator is a metal cylinder with an axial slot which may be 
used to adjust the inductance of a pulsed coil and to concentrate the magnetic field into a 
smaller volume. (See fig. 22.) The axial current components in the coil are antiparallel 
above and below the coil feed slot. The presence of the slot prevents these components 
from canceling each other out. The flux concentrator which is introduced into the coil 


32 




Figure 21.- Inductive lens for drift correction. 



Figure 22.- Diagram showing cross section of 
flux concentrator in coil. 

has its slot located on the opposite side from the coil feed slot. Thus the axial current 
components produced on the outer skin of the flux concentrator can cancel each other in 
closed loops. Only the azimuthal current components can reach the interior of the flux 
concentrator and be seen by the plasma. (See fig. 23.) This method not only reduces 
field nonuniformity due to current funneling but also, provided the radial thickness of the 
concentrator is not too large, guards against the cur rent- transfer effect. (The radial 
thickness of the flux concentrator should be large compared with the skin depth, however.) 
In fact, the flux concentrator should cancel any axial current which is symmetric about 
the median plane. Currents which cross the median plane may be caused by bank jitter; 
however, the sort of switch commonly employed has jitter of less than 50 nanoseconds, 
so these effects are presumably small. 
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Flux concentrator 



Figure 23.- Flux concentrator for shielding axial coil currents. 


THE FIELD-SHAPING EXPERIMENT 


A 1/2-scale replica of the Langley theta pinch has been built in order to measure 
the vacuum -magnetic-field distribution in the theta-pinch coil. Such an auxiliary facility 
is useful for evaluating the effect on field uniformity of various coil-collector geometries. 
In particular, the relative merits of collector-plate extension and flux concentrator can 
be compared. 


Description of Experimental Apparatus 

In order that the scaling include not only the geometry but also the skin depth, vari- 
able inductors are added in series with the collector plates so that the oscillation fre- 
quency of the coil current can be adjusted to approximately four times that of the full- 
scale experiment. Since the skin depth is proportional to w“V2 ) an increase of a factor 
of four in the frequency scales the skin depth down by a factor of two. 

The collector plates of the Langley magnetic compression experiment have been 
faithfully scaled in all details including the bolt holes and the slots for limiting short- 
circuited current. The presence of the lead bricks above the plates is simulated by a 
sheet-aluminum structure having the same shape as the lead bricks on the full-scale 
experiment. The lower lead bricks are simulated by aluminum plates bolted to the sup- 
port structure. The coils and flux concentrators are constructed of aluminum, which has 
electrical conductivity comparable to that of the beryllium -copper theta-pinch coil. Fig- 
ure 24 shows this field-shaping-experiment apparatus with a 70-centimeter-long coil 
made up of nine sections, each with a 5.5-centimeter bore. This figure should be com- 
pared with a similar view of the full-scale device shown in figure 10. Figure 25 shows 
the apparatus with the upper lead-brick simulator removed. The three variable inductors 
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Figure 24.- The l/2-scaIe model of the Langley theta pinch used for the field-shaping experiment. L-67-2086.1 



Figure 25.- Field-shaping-experiment apparatus with upper lead-brick simulator removed. L-67-2087.1 
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can be seen attached vertically to the back of the plates. Each inductor consists of a 
rectangular loop whose length may be varied. They are mounted to the collector plates 
by aluminum connector plates. 



The power source for the field-shaping experiment is a 37. 5-kilojoule capacitor 
bank composed of 30 1-microfarad clamshell capacitors. The capacitors are switched 
by modified commercial pressurized spark-gap switches. The trigger pulse for these 
switches is supplied by a 0.75-microfarad capacitor charged to 20 kilovolts. This capac- 
itor in turn is triggered by a resistor-capacitor (RC) circuit which is switched by a thy- 
ratron through a 5:1 pulse transformer to give a 35-kilovolt output. The measured volt- 
age rate of rise at the switch trigger pin as measured by a fast cathode ray oscilloscope 
(CRO) is greater than 1 MV/ps. This rise time insures a switch jitter of less than 
20 nanoseconds. The capacitor bank is shown in figure 26. It is located on a mezzanine 
directly over the experiment. The current is carried by one cable per capacitor over 


Figure 26.- Capacitor bank for the field-shaping experiment. 


L-67-2379 
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slack spools to the inductors. All cables are 5.1 meters long and have magnetically 
swaged terminations to reduce corona and arcing. 

Initiation of the firing signal by the control system triggers the pulse generator and 
the CRG camera shutters which are set for a 1 -second exposure. Approximately 0.5 sec- 
ond later the pulse generator triggers the bank and the oscilloscopes. 

The voltage on the capacitor bank is monitored by an electrostatic voltmeter. By 
approaching the final charging voltage slowly, repeatability in bank voltage can be kept to 
within 0.1 percent. This accuracy has been verified by checking the repeatability of 
probe traces on a CRO with the use of a type-Z preamplifier to raise the scale height to 
50 centimeters. In this way, changes from shot to shot of 0.1 percent could be detected. 

Use of Magnetic Probes 

Magnetic probes were used to measure small departures from uniformity of the 
magnetic field in the 1/2-scale theta pinch. In order to measure the local transverse 
force which would be exerted on a plasma column, one must measure the reciprocal of 
the local radius of curvature of field lines (see eq. (6)) given by 

1 _ VB 
R Bz 

9B r 8B Z 

where VB = — - = - — since the field is curl free. Hence, a precise point-by-point 
8z 8r 

measurement of B z and B r would determine 1/R. Theoretically one could just mea- 
sure B z and determine its radial gradient; however, this process is more difficult than 
using measurements of B r since the changes in B r are comparable with B r itself. 
The problem in measuring B r is that one is attempting to measure a small field com- 
ponent in the presence of a large field component. This difficulty can be surmounted, 
however. A probe oriented to pick up B r will generate a voltage proportional to 
B r + oqB z , where aq is the discrimination coefficient indicating the amount of B z 
pickup. When the probe is flipped 180° about the Z-axis, the voltage output becomes pro- 
portional to -B r + o^Bz, where i s the discrimination coefficient. The difference in 
the two signals is 2B r + (ctq - o^Bz which is equal to 2B r if oq = a^. The problem 
then is to make a probe which will maintain the discrimination coefficient a constant 
during the flip. Also oB z should be as small as possible in order to maximize the 
accuracy. The B z pickup comes from probe misalinement and (to a larger extent) from 
pickup in the leads. Therefore a probe was designed which would: (1) minimize stray 
pickup in the leads; (2) assure rigid leads so that lead loop areas would not change in the 
flip; (3) guarantee repeatable and accurate 180° flips; (4) be highly precise in alinement; 
and (5) have high sensitivity. The five design criteria were achieved in the probe shown 
schematically in figure 27. Stray pickup is minimized by using an inter-8-weave cable 
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Figure 27.- Schematic diagram of magnetic probe. 


with braid shielding which starts at the probe coil and continues via a copper cup and a 
uhf plug into the RG-52/U cable which runs to a screen room. In this way complete 
coaxial shielding is provided from the probe coil to the screen room and no stray loops 
are exposed. Additional shielding is provided by running the coaxial cable from the probe 
holder to the conduit in copper tubing. The leads are fixed in the teflon probe holder 
which has a 1 -centimeter -square cross section. The probe fits in a plexiglass probe 
carriage which has a 1- by 4.5-centimeter slot all the way through it for the probe holder. 
The probe is positioned radially by 5-millimeter teflon shims. The axial location is pro- 
vided in 10-centimeter steps by moving the probe carriage. Within the 10-centimeter 
steps, 1-centimeter increments are marked on the probe holder. This gives a radial 
positioning accuracy of ±0.1 millimeter and an axial positioning accuracy of ±0.2 milli- 
meter. The square cross section and precision fit of the holder, carriage, and shims 
provide a precise and repeatable flip angle of 180°. The probe arrangement inserted in 
the theta-pinch coil is shown in figure 28. The B r probe and a similar B z probe are 
shown in figure 29. The probe coils themselves are wound on a 6 -millimeter plexiglass 
rod, with 22 turns of enameled copper wire, The coils have a ratio of length to outer 
diameter of 0.72 which, according to reference 59, guarantees that the measurement is a 
point measurement and not an average over the fields occupying the coil volume. The 
coils are coated with a thin layer of polystyrene and carefully alined on the probe holder 
and fixed in place with a jig. The coils are then soldered to the leads and encapsulated 
in epoxy. 

In order to test the accuracy of the probe measurements, a Helmholtz coil was con- 
structed (see fig. 30). The Helmholtz coil, in principle, will produce a magnetic field 
whose line curvatures can be precisely calculated. The problem is that deviation of the 
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Figure 30.- The Helmholtz coll. 


L- 67- 2377 


field from perfection by, say, one part in a hundred due to disturbances from the leads is 
hard to avoid. This situation was shown to exist when a measurement of Br at the 
geometrical center of the coil, where B r should ideally be zero, showed that a value 
for Br of 1.28 percent of B z existed. Effort was made to reduce this interference, 
without appreciable improvement. The basic problem seems to be that the Helmholtz 
coil has a single turn. A multiturn coil would suffer less from lead interference; how- 
ever, its construction for high-voltage operation would be rather elaborate. Consequently, 
measurements of B r were made in the imperfect Helmholtz coil along an axis 5 centi- 
meters off center, from the median plane (z = 0) outward to the plane of one coil 
(z = 12 cm). These measurements are shown in figure 31. The height of the error bars 
is determined by the signal-to-noise ratio of the oscillograms. A curve was fitted to the 
points and then differentiated so that the field-line curvatures could be determined. Mea- 
sured values of radius of curvature are compared with values calculated for an ideal 
Helmholtz coil in figure 32. The close fit of the points to a smooth curve and the detec- 
tion of the small effect of negative curvature close to the median plane were highly 
encouraging. 
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Figure 31.- Measurements of B r in Helmholtz coil at r ~ 5 centimeters. 


Measurements of B r were made in the following six 1/2-scale theta-pinch 
configurations: 

(1) Configuration with normal coil and with lead -brick simulators but without 

extension (see fig. 24); 

(2) Configuration (1) with upper lead-brick simulator removed (see fig. 25) ; 

(3) Configuration with flux-concentrator coil and with lead-brick simulators but 

without extension; 

(4) Configuration (3) with upper lead-brick simulator removed; 

(5) Configuration with normal coil, with extension, and with lead-brick simulators 

placed 15.24 centimeters above extension to top collector plate; and 

(6) Configuration with flux- concentrator coil, with extension, and with lead-brick 

simulators placed 15.24 centimeters above extension to top collector plate 
(see fig. 33). 
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. L-67-2375.1 

Figure 33.- The I/2-scale model of the Langley theta pinch with extension and flux concentrator. 


The extension has a width equal to the coil length (70 centimeters) and places the 
coil 40.5 centimeters forward. The extension to the top collector plate is laminated into 
nine sections, each having the same width as the coil sections, with a 1.5-millimeter gap 
between them. The current flow in the extension to the bottom plate, even though it is 
unlaminated, follows the flow pattern in the extension to the top plate because of the skin 
effect. Hence, the effective length-width ratio of the extension is 5.4. 

The flux-concentrator coil is shown partially assembled in figure 34. The flux 
concentrator is a bored-out version of the normal coil such that the bore of the flux con- 
centrator, which has a wall thickness of 6.35 millimeters, is the same as the bore of the 
normal coil (5.5 centimeters). The concentrator and coil are machined from solid 
aluminum. 

For experimental measurement of the radial magnetic field in these six configura- 
tions, it was found that the signal-to-noise ratio could be optimized by reducing the 
capacitance of the bank to 6 microfarads and operating with a charging voltage of 10 kilo- 
volts. In order to maintain the frequency as near 80 kilocycles as possible (for purposes 
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Figure 34.- The flux-concentrator coil. 


L- 67 -2088.1 


of scaling the skin depth) , the inductors had to be tuned for maximum inductance. The 
result was an operating frequency of 96.7 kilocycles. The voltage on the bank appeared 
almost entirely across the inductors; only 250 volts were applied to the coil and thus the 
insulation problem was greatly simplified. 

Measurements of B r were made 7.5 and 17.5 millimeters offcenter radially in 
each direction in the plane of the coil feed slot. The radial direction toward the collector 
plate is taken as negative; the direction away from the coil slot is taken as positive. 
Measurements are made at intervals of 3 centimeters along the coil axis starting with 
z = 0 at the median plane and progressing to z = 30 centimeters, which is 5 centimeters 
from the end of the coil. The signals are carried to the screen room by a RG-52/U coax- 
ial cable and terminated with 50 ohms. A type-Z preamplifier with a CRO was used and 
the oscillograms were recorded photographically by an electrically actuated camera. The 
preamplifier features a dc bias voltage which can be used to raise the scale height and 
allow high gain to be used for maximum reading accuracy. This bias voltage was 
employed to give the signals a scale height of 6 centimeters whenever possible. The 
results of these measurements are presented in the following section. 

RESULTS AND DISCUSSION 
Field- Shaping-Experiment Results 

The results of the measurements of B r in the field-shaping experiment are shown 
in figure 35. 
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simulator but without upper lead-brick simulator or extension; flux 
concentrator has a 1-centimeter slot. 




(e) Configuration with flux-concentrator coil and with lead-brick simulators but with- 
out extension; flux concentrator has a 0.25-millimeter slot {oriented as shown in 
fig. 22). 


Figure 35.- Continued. 


47 



it) Configuration with flux-concentrator coil and with lead-brick simulators but with- (g) Configuration with flux-concentrator coil and with lead-brick simulators 
out extension; flux concentrator hjas a vertically oriented 0.25- mil ii meter slot. but without extension; flux concentrator has a 0.25- millimeter slot alined 

with the coil feed slot. 
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In figure 35(a) the gradient in B r is seen to diminish as the measurement location 
is moved away from the coil feed slot. This is essentially the same result found in ref- 
erence 4 for a 1/5-scale model. Ripples are observed in the B r /B z curve, the peaks 
of which correspond to the vertical slots between coil sections. The source of these 
ripples becomes evident upon inspecting figure 35(b), which shows results for a configu- 
ration that differs from the configuration considered in figure 35(a) only in that the upper 
lead-brick simulator has been removed. The amplitude of the ripples has been some- 
what reduced although the character of the curves is the same. This effect is due to cur- 
rents set up in the collector -plate structure. 

Figures 35(c) and 35(d) show the B r field inside the flux concentrator with and 
without the presence of the upper lead-brick simulator, respectively. Again, although 
the character of the curves is clearly the same, the ripple amplitude is somewhat less 
when the upper simulator is not present. Note that the more irregular field is on the 
positive side of the coil center line and that the field nearest the collector plate is quite 
smooth. Most important, the insertion of the flux concentrator has almost eliminated 
net field curvature. The effect of the ripples on the plasma is impossible to predict 
exactly, although it may well be that alternate regions of opposite curvature tend to 
inhibit instabilities by shorting out the charge separations which drive them. 

On the assumption that the rather irregular field near the flux-concentrator slot 
was caused by the large width of the slot (1 centimeter), a new flux concentrator was 
machined with a 0.25-millimeter slot and the experiment was rerun. The results are 
shown in figure 35(e). The amplitude of B r is reduced by approximately one -third and 
the variation of Br with axial distance close to the flux- concentrator slot simply sug- 
gests that some flux is lost as a result of smooth bending of field lines into the slot. 

These effects were confirmed by rotating the narrow-slot flux concentrator 90° such that 
the slot was vertical. The data taken for this configuration are shown in figure 35(f). 

This figure shows a very uniform field with an average amplitude for B r of less than 
0.1 percent of B z and no net gradient. In order to confirm the original concept of the 
flux concentration, the slot was next alined with the coil feed slot and measurements were 
made. It was anticipated that the result would be just as though no flux concentrator 
were present except that the ripples seen in figure 35(a), which were shown to correlate 
with the gaps between coil sections, would be absent. Figure 35(g) shows that the data 
obtained with the slots alined conform to the expected result. 

Figures 35(h) and 35(1) show the results for a collector -plate extension. These data 
reveal a highly uniform magnetic field and support the prior supposition that the current 
funneling is a dominating effect. It is interesting to note that the current-transfer effect 
(see section "Coil-Collector Discontinuity"), which was masked in figures 35(a) to 35(e) 
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and in figure 35(g) by either ripple or current funneling, is now in evidence. In fig- 
ure 3 5 (i) , for which the wide-slot flux concentrator is used, the field nearest the con- 
centrator slot (r = 17.5 mm) has a positive gradient. These effects are small and some- 
what localized and would probably not limit the plasma containment time. 

As seen from equation (6), the transverse force which causes plasma drift is pro- 
portional to 1/R, where R is the radius of curvature of the field lines. An average 
radius of curvature for each of the four radial positions in the six configurations has been 
calculated by considering values of B r only at z = 0 and z = 30. It is likely that the 
average curvature is the most influential factor in the drift characteristics of the plasma 
column. These average radii of curvature have been scaled up to apply to the full-scale 
Langley theta pinch (i.e., coil length of 140 centimeters and coil bore of 11 centimeters). 
The reciprocal of the radius of curvature 1/R is plotted as a function of radial position 
in figures 36, 37, 38, and 39, which show, respectively, the extension effect, the flux- 
concentrator effect, the lead-brick effect, and the extended-flux-concentrator effect. The 
direction of the forces is indicated by arrows on the data points. The selection of the end 
point at 30 centimeters for use in determining the gradient may be justified somewhat by 
the fact that it is neither in line with the vertical slot nor in the extreme end region where 
the field has begun to flare. The choice still has a certain arbitrariness about it, how- 
ever, and this, as well as the averaging of 1/R over z, must be weighed in evaluating 
the curves in figures 36 to 39. 

First, the effect of the extension is considered. A comparison is made between the 
scaled-up results for the normal coil without the extension (both with and without the 
upper lead-brick simulator present) and the scaled-up results for the normal coil on the 
extension. Figure 36 shows that the forces exerted when the extension is used are 
smaller than when it isn't. Moreover, the transfer effect causes the drift forces that do 
exist when the extension is used to be always directed toward the center of the coil rather 
than toward the wall. Therefore the extension should certainly suffice to stabilize the 
plasma column against drift. 

The flux-concentrator effect is shown next in a comparison of the scaled-up results 
for the flux-concentrator and normal coils both with and without the upper lead-brick sim- 
ulator present. The arbitrariness of the method by which the curves in figure 37 were 
obtained may have resulted in somewhat inconclusive results. The fact that the flux- 
concentrator radial field has such a large ripple makes the average curvature somewhat 
more uncertain than in the normal-coil cases. It would seem, however, that the average 
curvature is less for the flux-concentrator coil than for the corresponding normal coil. 

Unfortunately, the results are also inconclusive in evaluating the effect of the lead- 
brick simulator on the average curvature. The lead-brick- simulator effect is shown in 
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figure 38. For the normal coil the removal of the simulator appears to give a slight 
reduction in the drift forces, whereas for the flux-concentrator coil the presence of the 
simulator seems beneficial. In general, the previously stated conclusion, that the main 
effect of the large conducting mass on one side of the coil seems to be the introduction 
of ripple into the field, must stand as the only result which can be observed. 




Figure 38.- The lead-brick effect, as shown by a comparison of the scaled-up 
results for the presence and absence of the upper lead-brick simulator for 
both the normal and flux-concentrator coils. 
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The effect of an extended flux concentrator on the curvature is shown in figure 39. 
The average forces seem to be about the same for the normal coil on the extension and 
the flux-concentrator coil on the extension. The previously mentioned effect of the cur- 
rent transfer can be seen on both curves. 



Figure 39.- The effect of the extended flux concentrator, as shown by a com- 
parison between the scaled-up results obtained with the normal coil on the 
extension and the sea led -up results obtained with the flux-concentrator 
coil on the extension. 


Stabilizing Effect of Flux Concentrator on Drift in Full-Scale Theta Pinch 

Preliminary tests have been performed on the Langley theta pinch to determine the 
effect of a flux concentrator on the plasma drift. The flux concentrator used consisted of 
two pieces, with the center slot in an eight-section coil being left uncovered. The con- 
centrator was rolled from a 1.6-millimeter-thick copper sheet, with a 25-millimeter slot 
being left between the two pieces, and was insulated by poly (ethylene terephthalatej plastic 
film from the coil. 

Drift data were obtained with the slot of the flux concentrator on the same side of 
the coil as the current-feed slot and also with the flux concentrator rotated 180°. The 
results are compared in figure 40. In addition to streak photography of the plasma 
column by an image-converter camera, the intensity of strong lines of Si III at 3090 A 
was monitored photoelectrically side-on by a monochromator. The sweep duration of the 
streak camera was 50 ps. The filling pressure of the hydrogen in the tube was 100 milli- 
torr and the bank was operated at a charging voltage of 10 kilovolts to give an energy of 
250 kilojoules. 
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L-68-10,079 

Figure 40.- Drift control by flux concentrator. (The points on the streak photographs 
where the plasma strikes the tube wall are indicated by arrows.) 

One sees from the streak photograph on the left-hand side of figure 40 that the 
plasma struck the wall shortly before maximum current. This correlation is easily made 
as a result of the distortion in the image-converter tube optics produced by the field. 

Note also that the Si III signal rises strongly after 8 /is. The right-hand side of the fig- 
ure is for the flux concentrator oriented to stabilize the drift. The filling pressure and 
bank energy are identical to those used for the flux concentrator having its slot alined 
with the coil slot. The Si HI signal remains small during the whole test and the streak 
photograph shows a reasonably well-centered plasma column with only a slight suggestion 
of interaction with the inside wall near the end of the half-cycle of coil current. Unfor- 
tunately these results did not prove to be indefinitely repeatable. The stabilizing effect 
seemed to decay as additional shots were taken. Some of these shots, taken with the flux 
concentrator in the stabilizing position, are shown in figure 41 in chronological order. 
Progressive contamination of the tube wall may have been responsible for the breakdown 
of the stabilizing effect. Further work done on the plasma purity problem since these 
tests were made seems to bear out this hypothesis. 

It has been established that a definite correlation exists between pump-out time 
(a good indicator of the general cleanliness of the system) and the ability to produce a 
quiescent reproducible plasma. At the time the tests were made, considerable 
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L-68- 10,080 

Figure 41.- Decay of the flux- 
concentrator stabilizing effect prob- 
ably due to progressive wall con- 
tamination, as shown by streak 
photographs of 50-ps duration. 
(These photographs are presented 
in chronological order, with the 
one taken earliest at the top, and 
represent the trend through 
20 shots; the points where the 
plasma strikes the tube wall are 
indicated by arrows.) 


contamination existed and it is probable that axial motion in the turbulent plasma may 
have given rise to the decay of the drift stabilization as the contamination increased from 
shot to shot. 


CONCLUSIONS 


Arguments have been advanced to suggest that the major cause of drift in a theta- 
pinch plasma is magnetic-field-line curvature due mainly to the funneling of current from 
a long collector plate into a short coil. The rotational instability also seems to be influ- 
enced by magnetic-field uniformity although the exact mechanism is no doubt complex. 
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Magnetic-probe measurements have been made to ascertain the magnitude of the effect 
of the large conducting collector -plate structure, a collector -plate extension, and a flux 
coiicentrator on the vacuum -field structure in the coil. The flux-concentrator technique 
investigated by the author has been shown to reduce the field-line curvature, and tests 
performed on the full-scale theta pinch have suggested that a simple flux concentrator is 
capable of stabilizing the drift. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September 23, 1968, 

129-02-01-02-23, 
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